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ABSTRACT 

Correlations between the accretion luminosity and emission line luminosities (Lace 
and Ijiine) of pre-main sequence (PMS) stars have been published for many different 
spectral lines, which are used to estimate accretion rates. Despite the origin of those 
correlations is unknown, this could be attributed to direct or indirect physical rela¬ 
tions between the emission line formation and the accretion mechanism. This work 
shows that all (near-UV/optical/near-IR) Ijacc-^Une correlations are the result of the 
fact that the accretion luminosity and the stellar luminosity (L*) are correlated, and 
are not necessarily related with the physical origin of the line. Synthetic and observa¬ 
tional data are used to illustrate how the ^acc-^Une correlations depend on the Lacc-L* 
relationship. We conclude that because PMS stars show the Lacc-L* correlation im¬ 
mediately implies that Lace also correlates with the luminosity of all emission lines, 
for which the Lacc-Liine correlations alone do not prove any physical connection with 
accretion but can only be used with practical purposes to roughly estimate accretion 
rates. When looking for correlations with possible physical meaning, we suggest that 
Lacc/L* and L/i„e/L* should be used instead of Lace and L^ne- Finally, the finding 
that Lace has a steeper dependence on L* for T-Tauri stars than for intermediate- 
mass Herbig Ae/Be stars is also discussed. That is explained from the magnetospheric 
accretion scenario and the different photospheric properties in the near-UV. 

Key words: Stars: pre-main sequence-Stars: variables: T Tauri, Herbig Ae/Be- 
Accretion, accretion disks-Line: formation-Methods: miscellaneous 


1 INTRODUCTION 

The disk-to-star accretion rate is one of the most im¬ 
portant parameters driving the evolution of pre-main se¬ 
quence (PMS) stars. However, it is difficult to directly 
measure the mass accretion rate, for which indirect em¬ 
pirical methods are necessary to estimate it. A widely 
used method exploits the fact that the accretion luminos¬ 
ity (Lace) correlates with the luminosity of various emis¬ 
sion lines (Liine). Despite the unknown origin of these 
correlations, they are being used to quickly estimate ac¬ 
cretion rates. The Lacc-Lii„e empirical correlations have 

* E-mail: I.Mendigutia@leeds.ac.uk 


been derived using samples of PMS stars by comparing 
their accretion luminosities, mostly obtained from the UV 
excess and line veiling, with the emission line luminos- 


itv fsee e.e. iMuzerolle et al.l 

19983: iHerczee & Hillenbrandl 

l2008l:lDahmll2008l:lFane et alJ 

20091: iRieliaco et al.l 20121. and 


references therein). Currently, dozens of near-UV - opti¬ 
cal - near-IR spectral lines have been found to correlate 
with Lace for classical T Tauri (TT) stars (for instance, 
the hydrogen Balmer and Paschen seri es, Hel, 01, NalD 
and Call transitions, Bry... etc; see e.g. I Alcala et al.ll2of3 . 
AL14 hereafter). The correlations of the accretion lumi¬ 
nosity with several of these lines have been extended both 
to the sub-stellar and the intermediate-mass Herbig Ae/Be 
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(HAeBe) reRimes (|M ahjMTtv_et_alJ|200^; Rigliaco et al.ll201ll : 


iDonehew fc Brittain! 2011 ; iMendieuti'a et all 2011 , 2013a). 


Apart from the observational effort involved to look 
for additional emission lines that could serve as accretion 
tracers, several investigations aim to provide physical links 
between some of the spectral transitions and the accre¬ 
tion process, which wonld explain the origin of the Lacc- 
Liine correlations. In a nutshell, ei ther the lines are direct ly 
tracing the accreting region (e.g. Muzerolle et al.l Il998al lbl: 


iKurosawa et al.ll200B : [Rigliaco et al.ll2015l l. or they trace the 

accretion indirectly, b y probing the accretion-powered out¬ 
flows and winds (e.g. |_HartigarL_eLalJ 199^; Edwajds_et_alJ 


I 2 OO 6 I : IKurosawa et al.ll201ll : IKurosawa fc Romanovall20l3h 
The correlation with forbidden lines like |OI ] (6300 A) ex¬ 
hibited by HAeBes llMendiguti'a et ^|2011^ is more diffi¬ 
cult to explain, as this line is not identified with accre¬ 
tion/winds b ut rather with the surface layers of the circum- 
stellar disks dAcke et al.ll200^ . A further challenge to the 
varions explanations of the origin of the Lacc-Lune correla¬ 
tions is that the variations in the accretion rate as measnred 
from the UV excess do not generally c orrelate with the ob¬ 
served changes in th e line luminosities (Nguyen et al.| [2009l : 
ICostigan et al.lI 2 OI 2 I : iMendigutfa et ahl 2011 , 2013al l . How¬ 
ever, time delays between diffe rent physical processes could 
be present llDupree et al.l[^12ll . 


On the other hand, the accretion Inminosity is also 
found to correlate with the luminosity of the central star 
(L«). The Lacc-L* correlation extends over 10 orders 
of magnitnde in Lace, and ~ 7 orders of magnitnde in 


stellar to the HAeBe regime (see e.g. 

Natta et alJ 200GI: 

Clarke & Pringlell2006l: Tilling et al.ll2008 

: iMendigutfa et al.1 

2 OIII: IFairlamb et al.ll2015l. and references therein). Based 


snggested that the correlation between the accretion lumi¬ 
nosity and the luminosity of several emission lines in HAeBe 
stars could be driven by the common dependence of both lu¬ 
minosities on the stellar Inminosity. 


The main goal of this paper is to demonstrate the equiv¬ 
alence of the Lacc-Liine and Lacc-L* correlations. In partic¬ 
ular, we aim to show that all (near-UV, optical and near- 
IR) Lacc-Liine Correlations in PMS stars are driven by the 
relationship between the stellar Inminosity and the accre¬ 
tion luminosity, and that therefore the accretion luminos¬ 
ity necessarily correlates with the Inminosity of all spectral 
lines regardless of their physical origin. Section [2] introduces 
and partially re-analyses the Lacc-L» correlation in PMS 
stars. Section [3] shows the expression that links the Lacc- 
L« relationship with the Lacc-Liine correlations. The inter¬ 
dependence between both types of correlations is illnstrated 
in section[4]using both synthetic data and observational data 
from the literature. Some implications from all the previous 
analysis are included in section O Finally, section [6] snmma- 
rizes our main conclnsions. 


2 THE L/ 


7-l. correlation 


A representative example of the empirical correlation be¬ 
tween the accretion and stellar luminosities Q is shown in 
Fig. [ 1 ] It inclndes data from the literature for very low-mass 
TTs and sub-stellar objects/companions (log (Lt/L©) < - 
1.25), TTs (-1.25 < log (L«/L 0 ) < 0.75), late-type HAeBes 
(0.75 < log (L«/Lq) < 2.25), and early type HAeBes (log 
(L«/Lq) > 2.25). The sources belong to different star form¬ 
ing regions. The graph shows that Lace increases with L*, 
with a relation st eeper for the TTs than for t he HAeBes. 

According to lClarke fc Prin^ 1 I 2 OO 6 I ') and iTilling et al.1 
(l2008tl . the npper bound of the Lacc-L* correlation (Lace 
L,) is the consequence of sample selection effects; the 
luminosity of most stars above that limit is dominated 
by accretion and these objects are in a yonnger, em¬ 
bedded phase without an optically visible photosphere. 
The lower bound (Lace ~ O.OIL*, mainly for objects 
with L, > Lq) is limited by accretion detection thresh¬ 
olds (symbols with vertical bars in Fig. [T]). The physi¬ 
cal origin of the Lacc-L* correlation is the snbject of ac¬ 
tive debate. This topic is not analysed he re but we re¬ 
fer the reader to several related works fe.g. iPadoan et al.1 
2 OO 5 I: Alexander fc Armitagd 20061: Dullemond et al. 20061 : 


Vorobyov fc BasrJl2008l : lErcolano et al.ll2014h . Inst^d. our 

contribution below deals with the observed change in the 
slope of the Lgee-L* correlation between the TT and th e 
HAeBe stars (IMendigutfa et aI]|201ll : lFairlamb et al.l[2015l L 
We constructed a sample of artificial stars represent¬ 
ing the TT and HAeBe re gime by usin g synthetic mod¬ 
els of stellar atmospheres ( Kurued Il993tl . The properties 


of each object are provided in Table [T] Columns two and 
three show the stellar luminosity and effective temperature. 
From these, the stellar radii was derived, spanning between 
0.7 and 4 R© (column 4). The stellar masses (column 5) 
were derived assuming log g = 4, and cover the 0.2 - 6 M© 
range. Magnetospheric accretion (MA) shock modelling was 
carried out for each star by adding (blackbody) accretion 
contributions t o the photospheric (Kuruez) spectra (s ee e.g. 
the reviews in ICalvet et al.1 [ 2 OOOI : IMendigutfa l2013bl L Two 
representative examples are presented in Fig. [2] (left panel). 
The shock model was applied following the usual recipes 
fo r both the TTs and HAeBes, and we refer the re ader 
tolCalve t fc Gul lbririd (Il998l l ; IMendigutfa et al.] (l201ll ~) and 
iFairlamb et al. 1 2015h for further details. Three different val¬ 
ues for the UV excess in the Balmer reg ion of the spectra 
(from - 3500 to 4000 A, as defined in IMendigutfa et ahl 
l2013all were modelled for each object assuming typical val¬ 
ues for the inward flux of energy carried by the accretion 
columns (10^^ erg cm^ s“^) and the disk truncation radius 
(5R*): a ’’maximum” excess (0.70 magnitudes), whose cor¬ 
responding accretion contribution is Lace ~ L* for L« ^ 
Lq; a ’’minimum” excess (0.01 magnitudes) representative 
of the observational limit, and whose corresponding accre¬ 
tion contribution is Lace ~ O.OIL* for L* Lq; and finally, 
a ’’typical” excess in-between the two previous (0.12 magni¬ 
tudes). The resulting accretion luminosities are shown in the 
last three columns of Table [T] These are plotted versus the 


^ Its counterpart, the relationship between mass accretion rate 
and stellar mass, can be derived from the Lacc-L* correlation 


using PMS tracks (see e.g. lClarke &: Pringldl^OOSll . 
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log (L,/L ) 

® ^ ^ sun-^ 


Figure 1. Lacc— L* correlation for sub-stellar objects and TTs in 
differe nt star forming regions (crosses; with vertical bar s for upper 
limits: [Natta et ^l20od iHerczeg &: Hillenbrandll2008l , and refer¬ 
ences therein), four (s ub-) stellar/plane t ary compan i ons a round 
PMS stars (squares; IClose et alj 20141 : I Zhou et alJ l2014t) . the 
Lupus sample from AL14 (dark triang les), and HAeBes (cir¬ 
cles; with vertical b ars for upper limits: [Mendigutia et al]l2QlJ: 
iFairlamb et al.ll2015h . The red diagonal dotted lines indicate Lacc 
= L* and Lacc = O.OIL*. The three blue diagonal dashed lines 
represent the accretion luminosities expected from MA modelling 
for Balmer excesses of 0.70, 0.12 and 0.01 magnitudes (top, mid, 
and bottom lines, respectively). The vertical dotted line indicates 
the stellar luminosity at which the Balmer jump becomes appar¬ 
ent in the photospheric spectra (see also Fig. n. 


corresponding L* values (blue diagonal dashed lines in Fig. 
[!}, matching the overall distribution of data. We note that 
excesses larger than 0.70 magnitudes could still be measured 
for the less luminous sources (L* < Lq) without reaching the 
upper bound (Lacc ~ L*). 

The fact that the accretion luminosity increases with 
the stellar luminosity is a natural consequence of MA shock 
modelling. This is illustrated in the left panel of Fig. [51 When 
the same excess (flux ratio between the solid and dotted 
lines) is observed in stars of different stellar luminosity, the 
most luminous stars (blue dotted line) must necessarily have 
a larger accretion contribution (dot-dashed lines). In order 
to understand the different slope in the Lacc -L* correlation 
for TT and HAeBe stars, it is important to recall that the 
accretion contribution, and therefore Lacc, is proportional to 
both the temperature of the accretion columns (Tcoi) and 
the filling factor (/), which represents the fraction of the 
stellar surface covered by the accretion shocks. Variations 
in Tcoi and / move the accretion-generated continuum ex¬ 
cess along the wavelength axis and flux axis respectively. 
The typical value for Tcoi is ~ 10"^ K across the TT and 
HAeBe regimes. Therefore, the excess peaks close to the 
Balmer region for both types of star. However, their photo- 
spheric spectra (i.e. when accretion is not present) are signif- 


Table 1. Sample of artificial stars. Stellar parameters and accre¬ 
tion luminosities from MA. 


Star L* T* R,* M* (Lacc)m (Lacc)t (Lacc)M 

[log Lq] (K) (R©) (M©) [log L©] [log L©] [log L©] 


1 

-1.25 

3500 

0.65 

0.15 

-4.85 

-3.75 

-2.84 

2 

-1.00 

4000 

0.66 

0.16 

-4.28 

-3.17 

-2.26 

3 

-0.75 

4500 

0.70 

0.18 

-3.64 

-2.53 

-1.61 

4 

-0.50 

5000 

0.75 

0.20 

-3.10 

-1.99 

-1.04 

5 

-0.25 

5500 

0.83 

0.25 

-2.62 

-1.50 

-0.53 

6 

0.00 

6000 

0.93 

0.31 

-2.20 

-1.08 

-0.08 

7 

0.25 

6500 

1.05 

0.40 

-1.85 

-0.74 

0.29 

8 

0.50 

7000 

1.21 

0.53 

-1.57 

-0.45 

0.58 

9 

0.75 

7500 

1.41 

0.72 

-1.36 

-0.25 

0.76 

10 

1.00 

8000 

1.65 

0.99 

-1.13 

-0.02 

0.98 

11 

1.25 

8500 

1.95 

1.38 

-0.89 

0.22 

1.21 

12 

1.50 

9000 

2.32 

1.95 

-0.64 

0.47 

1.46 

13 

1.75 

9500 

2.78 

2.80 

-0.40 

0.72 

1.71 

14 

2.00 10000 

3.34 

4.05 

-0.15 

0.97 

1.96 

15 

2.25 

10500 

4.04 

5.93 

0.10 

1.22 

2.22 


Notes. Columns two to five show the stellar luminosity (loga¬ 
rithmic scale, from the integrated Kurucz model atmospheres), 
effective temperature, stellar radius and mass. Columns six to 
eight show the MA accretion luminosities (logarithmic scale) cor¬ 
responding to a minimum (m), typical (t) and maximum (M) 
Balmer excess of 0.01, 0.12 and 0.70 magnitudes, respectively. 


icantly different in that region. The Balmer jump becomes 
visible only for stars with log (L*/L©) ^ 0.25 (i.e. T* ^ 6500 
Kn This makes the spectra of stars with spectral types F 
and earlier more similar between them in the Balmer region 
than for later spectral types. Fig. [2] (top right panel) illus¬ 
trates the case; the photospheric U-B colour characterizing 
the Balmer region shows a steep dependence on the stellar 
temperature for cool stars, and flattens for hotter objects. 
Therefore, in order to reproduce a given Balmer excess, TTs 
require larger variations in the accretion luminosity than the 
ones that HAeBes need, for which the slope ALacc/AL* de¬ 
creases from the TT to the HAeBe regime. The accretion 
luminosity changes are mainly affected by variations in the 
filling factor. This is shown in Fig. [5] (bottom right panel), 
where the filling factors that are needed to reproduce the 
minimum, typical and maximum model excesses are plot¬ 
ted against the stellar temperature. The change of slope 
in this panel occurs at the temperature where the Balmer 
jump appears (-^ 6500 K), which corresponds to the stel¬ 
lar luminosity when the slope of the Lacc - L, correlation 
changes (log (L«/L 0 ) ~ 0.25). It is noted that we have ap¬ 
plied basic MA modelling without considering aspects like 
the chromospheric c ontribution to the spectra of TT stars 
llManara et al.ll2013l i or changes in the dis k truncation ra¬ 
dius depending on the stellar mass regime ( MuzgIoUg--£t_alJ 
I2OO4I : iMendigutfa et aP I2OIII : ICaulev fc .lohns-KrullI 2014l l . 
These factors could change the accretion estimates by less 
than 0.5 dex, without significantly affecting the modelled 
results in Figs. [T]and[ 2 ] 

In summary, the observed difference in the Lacc-L, cor¬ 
relation between TTs and HAeBes can be explained from 
the MA scenario and the differences in the near-UV stellar 


^ The Balmer jump disappears again in O stars with T* is 30000 
K 
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Figure 2. Left panel: MA modelling of a typical Balmer excess 
(0.12 magnitudes) for two representative stars with stellar tem¬ 
peratures of 7500 (blue) and 5500 K (red). The photospheric (Ku- 
rucz) spectra, the contribution from accretion, and the total flux 
obtained from the combination of the previous are represented by 
the dotted, dot-dashed, and solid lines, respectively. The fluxes 
are as they would be measured at the stell ar surface. Right pan¬ 
els: p hotospheric U - B colours (taken from iKenvon &: HartmannI 
Il995h characterizing the Balmer region of the spectrum (top) and 
Ailing factors necessary to reproduce a Balmer excess of 0.01, 0.12 
and 0.70 magnitudes (bottom) using MA, versus the stellar tem¬ 
perature. The vertical dotted line indicates the stellar tempera¬ 
ture at which the Balmer jump becomes apparent in the photo¬ 
spheric spectra (see also Fig.[J. 

properties between both types of stars. However, we em¬ 
phasize that the overall Lacc-L* correlation is not a mere 
consequence of the MA shock modelling but most probably 
reflects a deeper physical relationship between both parame¬ 
ters (see e.g. the references at the beginning of this section). 
For example, the specific slopes shown by different samples 
in different environments (see e.g the Lupus sample with 
solid triangles in Fig. [T]) cannot simply be explained from 
MA. Moreover, the Lacc-L* correlation seems to arise also 
in embedded, younger sources, when the accretion luminosi¬ 
ties are estimated from a variety of methods (Beltran & de 
Wit, to be submitted). Regardless of the underlying physical 
origin of the Lacc-L* correlation, for the rest of the paper it 
will be enough to remind that this arises whenever a signif¬ 
icant sample of PMS stars is considered. 


3 THE ACCRETION-STELLAR-LINE 
LUMINOSITY RELATION 

The relation between the accretion and stellar luminosities 
is usually expressed in the literature as Lace oc L*. This can 
also be written as a linear expression, which is a reasonable 
approach when the TT and HAeBe regimes are studied sep¬ 
arately. For a given star, we will assume that Lace and L* 
can then be related by: 

with a and b constants that depend on the star considered. 
When a sample of stars is studied, a and b represent the 
intercept and the slope of a linear fit to the data. This situ¬ 
ation will be analysed in section [4] 


The luminosity of a spectral line can be computed by 
multiplying the line equivalent width (EW) and the luminos¬ 
ity (per unit wavelength) of the adjacent continuum (L^): 

Liiac =LlxEW= X L., (2) 

with a the (dimensionless) excess of the (dereddened) con¬ 
tinuum with respect to the photosphere at the wavelength 
of the line (a = Ly/Z/A* ^ 1), and j3 the ratio between the 
total stellar luminosity and the stellar luminosity at that 
wavelength (/9 = L*/La* >> 1, in units of wavelength). The 
stellar luminosity in the second term of Eq.[2]was introduced 
in Eq. [T] obtaining: 

which is again a linear expression, with: 

A = a-6xlog(^^); 

B = b, 

where Lune/L* = aEW//?, is the line to stellar luminosity 
ratio. Therefore, if the accretion luminosity of a given star 
can be derived from its stellar luminosity through Eq. [T] 
then the same accretion luminosity can be recovered from 
the luminosity of any emission line through Eqs. [3] and U] 
with A and B constants that depend on the star and the line 
considered. Equations [T] and [3] are equivalent because both 
express a common dependence of the accretion luminosity 
on the stellar luminosity (Eq. [2}. 


4 THE DEPENDENCE OF THE Lacc-Bline 

CORRELATIONS ON THE Lacc-L* 

RELATION 

In this section we use both synthetic and empirical data to 
illustrate the dependence of the Lacc-Lune correlations on 
the Lacc-L* relation. Our first analysis provides a simple 
qualitative example on how the shape of the Lacc-L* rela¬ 
tionship has a strong effect on the Lace-Lane correlations. 
We use the sample of artificial stars introduced in the pre¬ 
vious section (see the first five columns of Table [T]). The 
Kuruez models were used to calculate L^ and /3 at 6000 A, 
whose values are presented in columns two and three of Ta- 
ble[2] Random EWs (between 1 and 10 A, column four) are 
assigned to each object. These range in EW is representative 
of emission lines with intermediate strength such as the Ca 
II or 01 transitions. The luminosity of an artificial emission 
line at 6000 A (column five) can then be obtained from Eqs. 

m 

The top left panel of Figure |3] shows two different Lacc- 
L* linear relations assumed for the sample. Both have the 
same intercept but a different slope. The reverse is shown 
in the bottom left panel, in which the slope is kept constant 
and the intercept varies. The right hand panels show the 
corresponding accretion luminosities versus the luminosity 
of the artificial line at 6000 A. The Lacc-Liine correlations 
follow the changes introduced in the Lacc-L* relation, vary¬ 
ing their slopes and intercepts. The range in the EW used 
only affects the scatter of the Lacc-Lune correlation, but this 
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Table 2. Sample of artificial stars. Continuum and line proper¬ 
ties. 


Star 

tt 

T ^ 

^6000 
[log L© A] 

/3(A = 6000 A) 

(A) 

EW 

(A) 

Leooo 
[log L©] 

1 

-5.64 

24317 

7 

-4.79 

2 

-5.13 

13363 

5 

-4.43 

3 

-4.74 

9790 

2 

-4.44 

4 

-4.42 

8343 

9 

-3.47 

5 

-4.14 

7792 

1 

-4.14 

6 

-3.88 

7573 

2 

-3.58 

7 

-3.63 

7587 

4 

-3.03 

8 

-3.39 

7698 

7 

-2.54 

9 

-3.15 

7880 

10 

-2.15 

10 

-2.92 

8284 

3 

-2.44 

11 

-2.69 

8773 

6 

-1.91 

12 

-2.48 

9603 

8 

-1.58 

13 

-2.27 

10581 

4 

-1.67 

14 

-2.07 

11723 

1 

-2.07 

15 

-1.86 

12905 

6 

-1.08 


Notes. Columns two to five list the luminosity of the continuum 
at 6000 A (logarithmic scale), the ratio between the total, star 
-h accretion, luminosity and the stellar luminosity at 6000 A, a 
random EW of an hypothetical emission line assigned to each 
star (between 1 and 10 A), and its corresponding luminosity at 
6000 A (logarithmic scale). 

is ultimately determined by the Lacc-L* relation. As intro¬ 
duced in section [3l the contribution of the continuum to the 
line luminosity dominates over the EW, and both the con¬ 
tinuum and the accretion luminosities are correlated with 
the stellar luminosity. In order to illustrate this, the EW 
range was increased multiplying by 10 all the EWs ^ 5 A in 
column 4 of Tabled and keeping the rest unmodified. This 
range in EW is representative of a strong emission line such 
as Ha. The new line luminosities are plotted with crosses 
in the bottom-right panel of Fig. [S] showing that for wider 
(narrower) EW ranges, the scatter in the Lacc-Liine corre¬ 
lation increases (decreases), but the correlation remains. 

Before using real data from the literature to illustrate 
how the Lacc-Liine empirical correlations are driven by the 
Lacc-L* relation, the equations described in the previous sec¬ 
tion have to be slightly modified. There, the values of Lace 
were given by Eqs. [T]and[3l where a, h, A and B differ de¬ 
pending on the individual star and spectral line. In practice, 
the values for the slopes and intercepts of these equations 
are estimated using linear regression fitting, which provide 
unique a and b values for a given sample of stars, as well 
as unique A and B values for a given spectral line. In this 
case it can be shown (see Appendix that the slopes and 
intercepts of the Lacc-L* and the Lacc-Lune empirical cor¬ 
relations are related by: 


B = bX€- (5) 

Inline X O’* 

e = - ~ 1, 

r* X Olijie 

where A, a; B, and b represent the intercepts and slopes 
of the Lacc-L* and Lacc-Liine correlations, as derived from 
least squares linear regression fitting, <log Liine/L*> the 



log (L„/L ) log (L /L ) 

® ^ * sun^ ® ^ 6000 sun^ 


Figure 3. Results for the sample of artificial stars, showing how 
changes in the assumed Lace—L* correlation (left panels) have an 
effect on the Lace—L^^^e relation (right panels) by changing the 
slope (top panels) and the intercept (bottom panels) of the former 
correlations. The crosses in the bottom-right panel represent the 
line luminosities obtained from a wider EW range, when the EWs 
^ 5 A in column 5 of Table [2] are multiplied by a factor 10. 

mean (logarithmic) line to stellar luminosity ratio, r* and 
Tiine the correlation coefficients of the Lacc-L* and Lacc- 
L^ine linear fits, and a* and aune the standard deviations of 
the log (L*/L 0 ) and log (L/^ne/L©) values. 

In short, when the empirical Lacc-L* and Lacc-L^ne 
correlations are compared, Eqs. should be used instead of 
Eqs. m These are slightly modified by including the e pa¬ 
rameter, which accounts for the fact that the empirical cor¬ 
relations are in practice derived from (least-squares) linear 
fitting. 

We use the observational data in AL14 to illustrate 
the dependence of the Lacc-Lii„e empirical correlations on 
the Lacc-L* relation. These authors studied a sample of 36 
low-mass TTs in the Lupus star forming region, for which 
they derived stellar parameters, accretion rates from the UV 
excess, and Lacc-Liine empirical correlations for dozens of 
emission lines in the spectral range from the near-UV to 

® Linear regression fits obtained from methods different than the 
usual least-squares are not considered in this work. The e pa¬ 
rameter should be eventually modified if other linear regression 
methods are used. 
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the near-infrared. To our knowledge, this work contains the 
largest number of spectral lines for which this type of corre¬ 
lations are derived. Another advantage is that for each star 
the accretion luminosity and the luminosity of all spectral 
lines were derived from the same spectrum, avoiding the 
problem of variability. In addition, all the stars are located 
at a similar distance, which guarantees that the correlations 
were not artihcially stretched when the fluxes are multiplied 
by the squared distances to derive the (accretion and line) 
luminosities. Therefore, we consider the Lacc-L, and Lacc- 
Liine correlations in AL14 as representative for similar cor¬ 
relations provided in the literature (see e.g the references in 
section [T|). 

The top panel of Fig. |4] shows the accretion and stel¬ 
lar luminosities of the stars studied by AL14. The observed 
trend is best htted by log (Lacc/L©) ~ -1.3 -I- 1.4xlog 
(L*/Lq) (solid line). The slopes and intercepts of the Lacc- 
hiine empirical correlations derived by AL14 (see their ta¬ 
ble 4), which are exactly recovered by Eqs. [S] are plotted 
in the mid and bottom panels of Fig. |4] versus e and e x 
<log Liine/L*>, respectively. The mid panel shows that the 
slopes of the Lacc-Lhne empirical correlations are a factor e 
smaller than the slope of the Lacc-L* correlation shown by 
the sample. As expected from Eqs. [S] the Lacc-Lune empiri¬ 
cal correlations become steeper when e increases, eventually 
reaching a slope of ~ 1.4 for e = 1. The bottom panel shows 
the expected linear decrease of the intercepts of the Lacc- 
Liine correlations with the (e-modified) line to stellar lumi¬ 
nosity ratio. Equations [5] also imply that the typical (mean) 
slope of all Lacc-Liine correlations is given by the slope of 
the Lacc-L* correlation of the sample, corrected by the mean 
value of e; <B> = 6 x <e>. Similarly, it can be derived that 
the mean intercept of the Lacc-Lune correlations is given by 
<A> = a - b X <£ X <log Liiae/L*>>. The two previous 
relations are also observed in the AL14 data, the mean val¬ 
ues indicated with the dashed lines perpendicular to both 
axis in the mid and bottom panels of Fig [5] 

In summary, the analysis of both a sample of artifi¬ 
cial stars and representative empirical data shows that the 
Lacc-Liine Correlations are driven by the underlying Lacc-L* 
relation shown by the sample of stars under study. 


5 CONSEQUENCES 

The first consequence of the analysis in the previous sections 
is that the fact that PMS stars show the Lacc-L* correla¬ 
tion immediately implies that Lace also correlates with the 
luminosity of any (near-UV-optical-near-UR) emission line, 
regardless of the physical origin of the spectral transition. 
Indeed, it even correlates with the luminosity of a randomly 
general artificial emission line (right panels of Fig. [^. As 
mentioned earlier, the scatter of the Lacc-Liine correlations 
increases when the lines’ EWs exhibit a larger range. A sim¬ 
ilar effect occurs for stars with strong excess at short, UV, 
wavelengths and long, IR, wavelengths. For lines observed 
a these short and long wavelengths, the ratio oEW/P (i.e. 
the line to stellar luminosity ratio; Eq. 13 ) becomes signif¬ 
icant, which could make the Lacc-Lune correlations much 
more scattered or eventually disappear. 

For the other lines, the Lacc-Lune correlations are 
mainly determined by the Lacc-L* dependence shown by 



log (L,/L ) 

® ^ * sun^ 





Figure 4. Based on results in AL14. Top panel: Accretion versus 
stellar luminosity. The best linear fit is log (Lacc/L 0 ) ~ -1.3 -h 
1.4xlog (L*/L 0 ) (solid line). Mid and bottom panels: Slopes and 
intercepts of the Lacc—^Hne empirical correlations versus the e pa¬ 
rameter (mid-panel) and the e-modified mean (logarithmic) line 
to stellar luminosity ratio (bottom-panel). The dashed lines indi¬ 
cate the mean values for the x and y axis, related from the slope 
and intercept of the top panel correlation by: <y> = 1.4<x> 
(mid panel), and <y> = -1.3 - 1.4<x> (bottom panel). 
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the sample under analysis. The intercepts and slopes pro¬ 
vided in the literature for the Lacc-L* correlation (a and 
b in Eg. [p vary depending on the sample of stars consid¬ 
ered (iFairlamb et al.l 120151 7 and references therein). Based 
on those works, a conservative observational limit is -2.5 ^ 
a ^ 0, 0.8 ^6^2. Conseqnently (see Eqs. |3]and[Sl), the 
slopes of all Lacc-Liine empirical correlations shonld also 
range in between ~ 0.8 and 2, whereas the intercepts shonld 
all be > 0 and decrease as the mean line to stellar Inminos- 
ity ratio increases. These predictions agree with all Lacc- 
Liine pnblished correlations based on observational data, to 
our knowledge. Interestingly, if two samples of stars show a 
different slope in their corresponding Lacc-L* correlations, 
then the slopes of the Lacc-Liine ones are simply related 
via B’ ~ B X (b’/b) (assuming that the e factors in Eq. [5] 
are roughly si milar in both samp l es). T his effect has already 
been observed. iMendieutfa et aP |20l3) reported a slight de¬ 
crease in the slope of the Lacc-L* correlation of a sample of 
34 HAeBe stars wit h respect to TTs (see also Fig. [1] and 
IFairlamb et al.|[2015l l. As discussed there, the slopes of the 
Lacc-Liine empirical correlations for the three lines studied 
(Ha, [01] (6300 A), and Bry) also show a similar decrease. 

That Lace correlates with hun^ is ultimately due to 
a common dependence of both luminosities on the stellar 
brightness. Because of this and the reasons above, the Lacc- 
Liine correlations alone cannot be seen as proof for either a 
direct or indirect physical connection between the spectral 
transitions and the accretion process. However, they are still 
useful expressions that can be applied to easily derive accre¬ 
tion luminosities without the need for sophisticated mod¬ 
elling of the UV excess. A basic measurement of a line lumi¬ 
nosity suffices. Given that both observational Lacc-Liine and 
Lacc-L* correlations show a roughly similar scatter (around 
± 1 dex in Lace), the latter can also be used to easily derive 
accretion rates from the stellar luminosity. 

Analogously, since Lii„e necessarily correlates with L* 
(Eq. ( 21 ), correlations between Lune and L* alone can not be 
taken as a possible physical link betwee n the spectral tran - 
sition and the stellar luminosity (see also lNatta et al.ll2014h . 
By extension, the luminosities of two different emission lines 
should also correlate with each other because of the common 
dependence on the stellar luminosity. Again, exceptions are 
possible for lines at sh ort/long wavelengt hs in stars with 
strong excesses (see e.g. iMeeus et al.l[2012l L 

In order to infer from correlations possible physical links 
involving the luminosity of a spectral line or the accretion 
luminosity, it is necessary to get rid of the common depen¬ 
dence of both parameters on the stellar luminosity. This can 
be done by dividing Lu^e and Lace by L*. Fig. [5] (top panels) 
shows the Lacc-Liine correlation for the sample of artificial 
stars from Table (21 and a given Lacc-L* relation, and the in¬ 
trinsic correlation between the stellar and line luminosities. 
However, the bottom left panels show that both Lacc/L* 
and L* do not correlate with Ljiae/L*, as expected from an 
artificial line created with random EWs. The right panels of 
the same figure show the results of the same exercise using 
real data from AL14. As expected, the Ha luminosity corre¬ 
lates with both the accretion and stellar luminosities, which 
as we have discussed has no possible physical interpretation. 
In contrast with the previous example, in this case the Ha 
line to stellar luminosity ratio is still correlated with the 
accretion to stellar luminosity ratio but not with the stellar 




log log (Lj^^/L.) 


Figure 5. Comparison between different luminosities normalized 
by the solar and the stellar luminosity, as indicated in the axes’ 
labels. The left panels refer to the sample of artificial stars in 
Table [J] and the right panels to real observations from AL14. 
Linear regression fits are overplotted for those cases with large 
enough correlation coefficients (r > 0.50; r-values indicated in 
each panel). 


luminosity itself, supporting the idea that this line is mainly 
driven by accretion and not by the stellar brightness. 

With this perspective in mind, we have confirmed that 
all line luminosities provided in AL14 correlate with each 
other, as expected. We also have checked that when the line 
luminosities are normalized by the stellar luminosities, some 
correlations remain while others disappear, indicating the 
presence or absence of a physical link between the different 
spectral transitions. For example, for Ha and Bry the cor¬ 
relation is not only between their line luminosities but also 
between their line to stellar luminosity ratios, suggesting a 
common physical origin for both transitions. In contrast, de¬ 
spite the fact that the luminosities of the Hell (4686 A) and 
the Call (8498 A) lines correlate, their line to stellar lumi¬ 
nosities do not show a significant correlation, suggesting a 
different physical origin. 

Finally, when the general Lace - L* correlation analysed 
in section (21 is transformed into Lacc/L* vs L*, no trend is 
shown either for the whole sample or for specific samples like 
the Lupus objects in AL14. The vast majority of the objects 
have 0.01 Bacc/L* ^ 1 (diagonal dotted lines in Fig. [TJ 
for all stellar luminosity bins. The typical value of Bacc/L* 
is 0.1, which corresponds to the modelled, typical Balmer 
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excess of 0.12 magnitudes. For the less luminous sources (L* 
< Lq), smaller Lacc/L* ratios can still be obtained from the 
same Balmer excess detection limit. As discussed in Sect.[2l 
this is the expected consequence of the MA scenario and the 
photospheric properties of the stars in the near-UV. 

It is beyond the scope of this work to carry out a de¬ 
tailed study on physical correlations involving stellar, line, 
and accretion luminosities. Instead, we have provided sev¬ 
eral examples to suggest that correlation analysis aiming to 
infer physical consequences should use Lii„e/L, and Lacc/L* 
and not simply Luns and Lace- 


6 SUMMARY AND CONCLUSIONS 

The Lacc-L* empirical correlation in PMS stars has been 
partially re-analysed taking into account the newly avail¬ 
able accretion rates for HAeBes. Despite the physical origin 
of the Lacc-L* correlation remains subject to debate, the 
observed change of slope from the TT to the HAeBe regime 
can be understood from the MA scenario and the near-UV 
photospheric properties of the stars. 

We have shown that the fact that PMS stars show the 
Lacc-L* correlation immediately implies that Lace also cor¬ 
relates with the luminosity of any (near-UV, optical, near- 
IR) emission line, regardless of the physical origin of the 
spectral transition. The overall Lacc-Lune trends are mainly 
governed by the Lacc-L* correlation shown by the sample of 
stars under analysis. In particular, the slopes of the Lacc- 
Liine empirical correlations should typically be between ~ 
0.8 and 2 for all spectral lines, which are the observational 
limits for the slope of the Lacc-L* relation. The intercepts 
also depend on the Lacc-L* correlation, all of which are > 
0 and increasing as the line to stellar luminosity ratio de¬ 
creases. 

Despite the fact that the Lacc-Lune correlations alone 
do not constitute an indication of any direct or indirect phys¬ 
ical link between the spectral transitions and accretion, they 
are a useful tool to easily derive estimates of the accretion 
rates. The Lacc-L* correlations can be used for the same 
purpose. Similarly, correlations between stellar and line lu¬ 
minosities, or between different line luminosities, do not in¬ 
dicate a physical relation between the parameters involved. 
Instead, we suggest that the line to stellar and accretion to 
stellar luminosity ratios should be used when investigating 
the possible physical origin of the various correlations. 
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of the hacc - L* and Lace - Liine linear correlations: 


B — t X b', 

f’line X O'* (A5) 

e = -. 

r* X Cline 

On the other hand, the mean value <log Lacc/LQ> can 
be found in the expression for a of Eq. IA2I and introduced 
in the expression for A of Eq. IA4I Also considering Eq. IA51 
the expression that relates both intercepts is: 


A = a—bxex 




X 



(A6) 

The third term could been neglected ((1 - e)/e ~ 0) com¬ 
pared with the two other terms in the previous equation. 


APPENDIX A: RELATION BETWEEN THE 
Lacc—Lline AND Lacc-L* linear 
REGRESSION CORRELATIONS 

Consider a sample of N stars for which measurements of 
accretion and stellar luminosities [log (Lacc/L 0 )i,..., log 
(Lacc/L 0 )iv; log (L*/L 0 )i,..., log (L*/L 0 )iv)] are available. 
A linear fit to the data provides an expression that links 
both variables through 

with a and b constants representing the intercept and the 
slope, which from least-squares linear regression are given 
by 



where r* is the correlation coefficient 1 for well correlated 
data), and aacc, o-»; <log Lacc/LQ>, and <log L«/L 0 > the 
standard deviations and the means of the log (Lacc/L 0 )i 
and log (L*/L 0 )i values, respectively. 

Similarly, if for the same sample of stars there are ad¬ 
ditional measurements of the luminosity of a given emission 
line [log (Liine/L 0 )i,..., log (Lii„e/L 0 )Ar], then a linear fit 
provides 

log(T^) =^ + ^^log(^) ’ (^3) 

with A and B constants given by least-squares linear regres¬ 
sion 



where the correlation coefficient, standard deviations, and 
means now refer to the [log (Lacc/L 0 )i, log (Liine/L 0 )i] val¬ 
ues. 

The standard deviation Uacc can be found in the expres¬ 
sion for b of Eq. IA21 and then introduced in the expression 
for B of Eq. IA4I providing the expression relating the slopes 










